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A DIFFERENTIAL CORRECTION TECHNIQUE FOR DETERMINING 
PHYSICAL PROPERTIES OF THE MARTIAN ATMOSPHERE BY USE OF 
SOLAR OCCULTATION AS SEEN FROM AN ARTIFICIAL SATELLITE 

By James R. Williams and H. Andrew Wallio 
Langley Research Center 

SUMMARY 

A satellite in a circular orbit about Mars is assumed to be equipped with instru- 
mentation capable of observing the diminution of sunlight passing through the atmosphere 
at grazing incidence to the Martian surface. A differential correction technique utilizing 
a least-squares method of data fitting has been used in conjunction with a model of the 
Martian atmosphere to determine the surface pressure, surface temperature, and tem- 
perature gradient. In the visible spectrum a wavelength of 5500 A (1 A = 10-10 m) pro- 
duces the quickest convergence for the widest range of varied initial conditions of the 
three parameters under consideration. 

INTRODUCTION 

Before a scientific payload can be soft-landed on Mars at precise locations, a 
better knowledge of the physical properties of the Martian atmosphere is required. One 
technique which could serve as an independent method in determining such atmospheric 
properties as pressure scale height, surface pressure, surface temperature, and tem- 
perature lapse rate (defined as the negative of the temperature gradient) is outlined in 
reference 1. This technique utilizes a photometer aboard an artificial satellite to mea- 
sure the variation in light intensity from the sun as the light passes at grazing incidence 
through the Martian atmosphere. 

Reference 1 demonstrates that the occultation technique allows ready discernment 
of differences in the curves of light intensity variation for various models of the atmo- 
sphere. The present paper investigates the problem of extracting the atmospheric prop- 
erties of surface pressure, surface temperature, and temperature lapse rate if a mea- 
sured curve of the light intensity variation is given. The purpose of this study is to 
determine whether a first-order differential correction technique is adequate to deter- 
mine these physical properties of the Martian atmosphere. 



In order to accomplish this a Taylor expansion of the equation of light intensity 
variation is truncated after the first-degree terms, and a least-squares fit is made to 
the data on the basis of an initial estimate, to determine the increments for the atmo- 
spheric properties that are to be used as the input conditions for the second iteration. 
This process is repeated until convergence criteria are met. 

As in reference 1, the factors which determine the mathematical model for the 
light intensity variation are hydrostatic equilibrium in the atmosphere, refraction through 
the atmosphere, limb darkening and finite size of the sun, and Rayleigh scattering in the 
atmosphere. A model using these factors appears to be adequate since recent findings 
of Mariners VI and VII (ref. 2) indicate that clouds and fog are not prominent in the 
Martian atmosphere, and therefore multiple and Mie scattering can be neglected. 

The technique presented in this paper can be used as a method to determine sur- 
face pressure and temperature, pressure scale height, and temperature lapse rate in the 
atmosphere. The advantage of using this occultation data is that with an attitude control 
system and cruise orientation similar to that used on the Lunar Orbiter and Mariner 
spacecraft, one spacecraft axis would be nominally pointed toward the sun, and a prop- 
erly mounted solar sensor would require no additional attitude control. 

SYMBOLS 


A 

D 

E 

g 

h 

I(x) 

lo 


matrix of partials 

distance from satellite to planet, 8700 km 
error 

acceleration due to gravity, 3.7 m/sec2 
pressure scale height, TR/mg, km 

intensity of sunlight passing through Martian atmosphere 
intensity of sunlight in free space 


*x 

ratio of 

I(x) 

Ix ,i 

value of 

h 


2 



l( y) factor determining limb darkening and finite sun size 

m average mass of gas, kg/mol 

p surface pressure, mb 

R universal gas constant 

r radius of Mars, km 

T surface temperature, K 

V variable limit of integration 

x coordinate of satellite trajectory, km 

y dummy variable of integration 

z altitude above Martian surface, km 

z max maximum value of z 

a half-angle subtended by sun as seen from satellite, radians 

0 Rayleigh scattering coefficient, km'l 

r temperature gradient parameter, km'^ 

X wavelength, A (1 A = 10“ 10 m ) 

co deviation of light ray from straight line at Martian surface, radians 

Primed parameters are estimated values. Parameters prefixed with A represent 
the difference between data and estimate. 

ASSUMPTIONS AND ANALYSIS 

In this paper the major effects on the variation of the light intensity are differential 
refraction through the Martian atmosphere, limb darkening of the sun, and Rayleigh 


3 


scattering due to the gas molecules in the atmosphere. The mathematical model of these 
effects is derived in reference 1. Mie scattering and multiple scattering due to clouds 
or fog are neglected (see ref. 2). 

Differential refraction is the divergence of two infinitesimally separated parallel 
rays which traverse an atmosphere whose index of refraction is a function of altitude. 

Limb darkening of the sun occurs because the light which reaches the satellite is 
the sum of rays from the various parts of the disk of the sun, each ray passing through a 
different layer of the Martian atmosphere. The effect is a spreading of the curve of 
light intensity variation and a smoothing of the discontinuities due to the fine structure 
of the atmosphere. Implicit in limb darkening effects is the consideration of the finite 
size of the sun. 

Rayleigh scattering is the scattering of light by the gas molecules in the atmosphere, 
where the wavelength of the incident radiation is long in comparison to the dimensions of 
the molecule. The scattering is dependent on the fourth power of the wavelength and the 
number and type of gas molecules in the atmosphere. 

Thus, for reference 1, the ratio of the light intensity of the sun received through 
the atmosphere to the light intensity of the sun in free space is given by combining the 
effects of refraction, limb darkening, and Rayleigh scattering as follows: 


Io 




with the auxiliary equation 


(° <z<z max) (1) 


x + aDy = -z + wD ^1 + ^ Th - Tzjexp^^ 
and the variable limits of integration given by 


(2) 


fi 


V = 


ckD 


a>D 


1 + ? 



X 


l - 1 


( la HsR 1 + i rh H) 

(- 1 -^[“ D ( 1 + I rh )- X ] sl ) | 

(- i5 ^H i+ f rh H) j 


( 3 ) 
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For brevity of notation, the following definition is introduced: 


1* = ^ < 4 > 

A o 

In this paper it is assumed that the atmosphere is composed of 90 percent carbon dioxide 
and 10 percent argon. (See ref. 2.) The remaining assumptions are the same as in 
reference 1. 

The dashed curve in figure 1 represents a typical variation of light intensity with 
distance into the geometric shadow of the planet. Since the constituents of the Martian 
atmosphere have been assumed, I x is a function of r, p, and T only: 


Ix = yr^T) 


(5) 


The functions r, p, and T are the unknown physical parameters that characterize the 
Martian atmosphere. The technique used in this paper to solve for these parameters 
consists of assuming a light curve which approximates the observed light curve and for 
which the parameters r f , p f , and T f are known. (Primed parameters are estimated 
values.) Now from Taylor’s theorem, 


I x ( r >P>T) 


V (r’, P ’,T’) + 


9l x , 

^L(r 


r’) +^-( P - p') +-%-(t - t’) + . . . 

9 p, vr r 9T , 


( 6 ) 


where second and higher order terms have been dropped. By definition, 


> 

1 

II 

< 

(7) 

at = r - r* 

(8) 

Ap = p - p’ 

(9) 

AT = T - T' 

(10) 


Equation (6) may be rewritten in a matrix notation as 

/ ar \ 

Ap 


AI, 


9l x' aI x' 9l X* ' 

ar r 8p* 8T' / 


W / 


( 11 ) 
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For brevity let 


A = 


ay s y #y\ 

ar* ap aT* / 


( 12 ) 


Solving equation (11) |under the constraint that 
and AT yields the least-squares solution 


2X* 


is a minimum 


for AT, Ap, 



The solid curve in figure 1 is a typical first estimate y of I x , the dashed 
curve. Now A^ is computed at approximately 100 points from figure 1. These values, 
along with the elements of A, are substituted into equation (12) and values for AT, 

Ap, and AT are computed in a least-squares sense. These results are substituted into 
equations (7), (8), (9), and (10) to obtain new values for T, p, and T which are used 
to generate a new estimate for the I x ' curve. The process is repeated for either 
41 iterations or until the error is less than 10"^. The error is defined as 

100 

E - 1 K,i) 2 

i=l 

where AI X ^ is computed at the ith point. 

DISCUSSION AND RESULTS 


The values chosen for the assumed measured light-intensity curve are T = 240 K, 
r = -5 X 10 km - •*■, and p = 10 mb (1 mb = 100 N/m2). The surface temperature was 
arbitrarily chosen between that of an illuminated hemisphere (252 K) and a planet average 
theoretical temperature (213 K) calculated for solar radiant energy equilibrium (ref. 3). 
The surface pressure was chosen to be between a low of 5 mb measured by Mariner IV 
(ref. 2) and the value of 17 mb calculated from spectroscopic observation of the 
1.6 pm CO2 bands (ref. 4). The value of r, which is related to the temperature lapse 
rate -TT, was chosen between the values of -3 X 10 - ^ km - * and -6.42 X 10"^ km“l 
obtained from a least-squares fit to the temperature curves of the Martian atmospheric 
engineering models used by the Viking Project Office. 
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For the initial conditions, only one of the three parameters was varied while the 
other two were started at the desired convergence values. The question of varying the 
three parameters simultaneously for input is essentially considered, because in the sec- 
ond iteration all three parameters are displaced from the desired conditions. Therefore, 
if the second iteration is considered as the first, all three parameters are varied 
simultaneously. 

Reference 1 pointed out that the shorter visible wavelengths (3500 A) produced the 
most readily discernible effects of the light intensity variation. Because of the impor- 
tance of wavelength and its effect on the results, this paper considers four wavelengths: 
3500 A, 4700 A, 5500 A, and 6500 A. 

The results are presented in tables and graphs. The tables list the final conver- 
gence values of p, T, and T along with their relative errors, and thus provide a quick 
summary of the results. On the other hand, the graphs of p’, T', and r* plotted 

against iteration number provide insight into the rate of convergence to the correct values 
and the convergence to local minimums. 

In figures 2 to 6 the initial p' values range from 2 to 25 mb while the other initial 
values are T' = 240 K and F = -5 x 10"^ km“l. The parameters are plotted against 
iteration number for the four wavelengths under consideration. 

Table I summarizes the results of figures 2 to 6. 

Figures 2 and 3 illustrate the rate of convergence when p is estimated lower than 
the desired value. Figure 2 is for an initial value of p' = 2 mb, and figure 3 is for an 
initial value of p’ = 6 mb. Figures 2(a), 2(b), and 2(c) converge within 3 percent of the 
desired values of p and T after 40 iterations. (See table I.) The desired value of T 
is being approached in figures 2(a), 2(b), and 2(c); however, its convergence rate is much 
slower than that of p and T. The temperature gradient parameter T is slowly con- 
verging in figures 2(a) and 2(b) but after 40 iterations its value is too great, by 20 percent 
and 25 percent, respectively. (See table I.) In figure 2(c), for X = 5500 A, the value 
of r after 40 iterations is within 1 percent. (See table I.) The same statements as 
were made about figures 2(a), 2(b), and 2(c) may be made about figures 3(a), 3(b), and 3(c), 
with 3(c) showing the best convergence after 40 iterations from the point of view of atmo- 
spheric physics. Figures 2(d) and 3(d), for X = 6500 A, do not converge to physically 
acceptable values of T and T. 

Figures 4, 5, and 6 are plots of the convergence rate for initial p’ values of 
14 mb, 18 mb, and 25 mb, respectively. Except in figure 6(d), all the values are con- 
verging for initial estimates of p' which are too high; however, the case for 
X = 4700 A converges most slowly, as can be seen in figures 4(b) and 5(b). In figure 6(d) 
the three parameters are far from correct values after 40 iterations. For the three 
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cases studied with p T initially estimated too high, the wavelengths that converge best 
are 5500 A and 3500 A. (See table I.) These results, together with the results for the 
cases in which p' was estimated too low, indicate that 5500 A is the most acceptable 
wavelength for determining an accurate p by using a linear estimation technique. 

In figures 7 to 10 the initial T’ values range from 160 to 320 K while the other 
initial values are p’ = 10 mb and T’ = -5 X 10"^ km -1 . The parameters are plotted 
against iteration number for the four wavelengths under consideration. 

Table II summarizes the results of figures 7 to 10. 

As was true of the first set of figures, T is the parameter which is the most dif- 
ficult to converge. Figures 7 and 8 are convergence plots for the cases of T’ esti- 
mated lower than T; here T’ = 160 K and 200 K, respectively. For a T' of 160 K the 

o 

best convergence for r is at a wavelength of 5500 A. At the end of 40 iterations, fig- 
ure 7(c) has r converged within 11 percent, p within 4 percent, and T within 8 per- 
cent. For these percentages and the final values of p, T, and T, see table II. For 
T' = 160 K, the p and T values after 40 iterations are better at 3500 A than at 5500 A, 
but r converges only within 32 percent at 3500 A. (See table II.) In figure 8, for 
T' = 200 K, there is convergence within 2 percent of the desired value for all three 
parameters at both 5500 A and 6500 A. 

Figures 9 and 10 are convergence plots for T' estimated higher than T; here 
T f is 280 K and 320 K, respectively. For a T' of 280 K there is convergence for all 
four wavelengths, with the best convergence for all three parameters occurring at 5500 A 
as shown in figure 9(c). Figures 10(a), 10(b), and 10(c) show convergence for all param- 
eters to within 10 percent of the desired value, but in figure 10(d) the parameters fail to 

o 

converge. (See table II.) However, a wavelength of 5500 A, shown in figure 10(c), again 
produces the quickest convergence for a given iteration. For a surface temperature that 

o 

can be estimated within 33 percent (see table II), a wavelength of 5500 A appears to be 
most suitable for obtaining convergence of p, T, and T. 

In figures 11 to 15 the initial T' values range from -50 x 10"^ km"l to 
5 X 10 km - * while the other initial values are p' = 10 mb and T' = 240 K. The 
parameters are plotted against iteration number for the four wavelengths under 
consideration. 

Table HI summarizes the results of figures 11 to 15. 

Figure 11 presents plots for T' = +5 X 10"3 km - *, which is of the desired magni- 
tude but has the wrong sign. Only the case with \ = 5500 A (fig. 11(c)) converges 
within 3.5 percent for all three parameters. (See table HI.) Figure 12 presents plots 
for r* = 0 as the initial value. For this case both X = 5500 A and \ = 6500 A con- 
verge within 2.6 percent for all three parameters. (See table HI.) Figure 13 is for an 
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initial T’ of -8 X 10 km - -*-. For this case all parameters converge at all wavelengths 

except X = 3500 A, and even at X = 3500 A only T, with an error of 11.2 percent, is 
not close to the desired value. (See table in.) Figure 14 is for T' = -12 X 10 - ^ km - * 
as the initial value. In figure 14 only the parameters at X = 3500 A do not converge to 
the desired values. In figure 15 the initial estimate of I” is in error by an order of 
magnitude, or T’ = -50 X 10 km - *. The only wavelength for which convergence is 
attained for all the parameters after 40 iterations is 5500 A. Thus, for an estimate 
of r within 900 percent a wavelength of 5500 A appears to be most suitable for conver- 
gence of p, T, and r. 

The results of the convergence figures indicate that a wavelength of about 5500 A 
should give the best overall convergence for a wide range of input conditions, when the 
linear estimation technique is used. At this wavelength the best convergence is obtained 
for p, the next best for T, and the worst for r. Thus, this technique is most sensitive 
to p, T, and T in that order. This is an interesting fact because in other investiga- 
tions these parameters have been best observed and measured in the inverse order. In a 
great many studies T is taken, to a first approximation, as being zero; that is, an iso- 
thermal atmosphere is assumed. The surface temperature can be estimated fairly well 
by assuming thermal equilibrium of the planet and the sun, and it has also been well mea- 
sured from polarization studies. The property that has yielded the widest range of 
measurements is the surface pressure, which varies from a calculated 5 mb (Mariner IV) 
to a calculated 87 mb (ref. 5) . 


CONCLUDING REMARKS 

A first-order differential correction technique utilizing a least-squares method of 
data fitting has been applied to the sun occultation light curves to recover the surface 
pressure, surface temperature, and temperature lapse rate in the Martian atmosphere. 
The model of light passage in the atmosphere included the effects of refraction, limb 
darkening, finite sun size, and Rayleigh scattering. The results of this investigation 
demonstrate that the most suitable wavelength at which to obtain data for quick conver- 
gence, even when the initial estimates are in error by 150 percent for surface pres- 
sure p, 33 percent for surface temperature T, or 900 percent for temperature gradient 

o 

parameter T, is in the 5500 A region of the visible spectrum. At this wavelength, the 
order of best convergence was found to be (1) p, (2) T, and (3) T. 
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This technique of utilizing the occultation of sunlight as seen by an artificial satel- 
lite and the differential correction technique to recover atmospheric parameters can be 
used as an independent method to determine the Martian atmospheric parameters or to 
provide a check on any other data source. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., March 29, 1971. 
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TABLE I.- SUMMARY OF RESULTS FOR VARIED INITIAL VALUES OF p' 


Figure 

Wavelength, 

A 


Initial conditions 


Final conditions 


Error, percentage 

,P', 

mb 

T' 

,K 

F, km-1 

p, mb 

T, K 

T, km' 1 

p 1 ‘ 10 v ioo 

T* - 240 v x 

r + (5 X io-3) v 

10 

240 

-5 X 10-3 

2(a) 

3500 


2 

240 

-5 X 10-3 

9.97 

240.4 

-5.98 X 10-3 

0.3 

0.2 

19.6 

(b) 

4700 







9.87 

243.6 

-6.27 X 10-3 

1.3 

1.5 

25.4 

(c) 

5500 







9.99 

246.8 

-5.05 X 10~ 3 

.1 

2.9 

1.0 

(d) 

6500 

\ 

/ 





31.1 

1634.6 

-4.55 X lO- 1 ^ 

211.0 

581.1 

100.0 

3(a) 

3500 


6 





9.97 

240.5 

-5.83 X 10-3 

.3 

.2 

16.6 

(b) 

4700 







9.90 

241.8 

-5.67 X 10-3 

1.0 

.8 

13.4 

(c) 

5500 







9.99 

242.4 

-5.07 X lO- 3 

.1 

1.0 

1.4 

(d) 

6500 

N 

/ 





31.1 

1634.6 

-4.55 X 10- 15 

211.0 

581.1 

100.0 

4(a) 

3500 

14 





10.0 

240.1 

-4.87 X lO" 3 

0 

.1 

2.6 

(b) 

4700 







10.1 

238.2 

-4.07 X 10" 3 

1.0 

.7 

18.6 

(c) 

5500 







10.1 

236.0 

-4.94 X 10-3 

1.0 

1.7 

1.2 

(d) 

6500 

\ 

/ 





9.99 

238.6 

-5.10 X lO' 3 

.1 

.6 

2.0 

5(a) 

3500 

18 





10.0 

240.0 

-4.44 X lO' 3 

0 

0 

11.2 

(b) 

4700 







10.2 

237.4 

-3.61 X lO" 3 

2.0 

1.1 

27.8 

(c) 

5500 







10.2 

231.5 

-4.81 X lO' 3 

2.0 

3.5 

3.8 

(d) 

6500 

\ 






9.98 

237.5 

-5.16 X lO" 3 

.2 

1.1 

3.2 

6(a) 

3500 

25 





10.0 

240.1 

-4.75 X lO -3 

0 

.1 

5.0 

(b) 

4700 







10.1 

239.0 

-4.48 X lO -3 

1.0 

.4 

10.4 

(c) 

5500 







10.2 

227.2 

-4.66 X lO -3 

2.0 

5.3 

6.8 

(d) 

6500 

> 

/ 

V 

/ 

\ 

/ 

10.9 

159.2 

-14.2 XIO'3 

9.0 

i 

33.7 

184.0 


TABLE H.- SUMMARY OF RESULTS FOR VARIED INITIAL VALUES OF T' 


Figure 

Wavelength, 

A 

Initial conditions 

Final conditions 

Error, percentage 

p f , 

mb 

T > 

K 

r’, km-l 

p, mb 

T, K 

T, km-1 

c 

c 

T" 

> 

O 

i— 1 

I 

T ’ - 240 v 1()0 

r + (5 x io-3) v 100 

10 

240 

-5 X lO" 3 

7(a) 

3500 

10 

160 

-5 X 10" 3 

10.0 

239.8 

-3.40 X 10 -3 

0 

0.1 

32.0 

(b) 

4700 







10.2 

237.4 

-3.58 x 10-3 

2.0 

1.1 

28.4 

(c) 

5500 







10.4 

221.2 

-4.47 X IQ' 3 

4.0 

7.8 

10.6 , 

(d) 

6500 



\ 

/ 



9.97 

196.4 |-9.13 X 10-3 

.3 

18.2 

82.6 

8(a) 

3500 



200 



10.0 

239.8 

-3.51 X 10-3 

0 

•1 

29.8 

(b) 

4700 







10.1 

237.6 

-3.77 x 10-3 

1.0 

1.0 

24.6 

(c) 

5500 







io.i 

236.5 

-4.92 X 10-3 

1.0 

1.5 

1.6 ; 

(d) 

6500 

1 

i 

' \ 

/ 



10.0 

238.5 

' -5.10 x 10-3 

i 

0 

.6 

2.0 

9(a) 

3500 



280 



i 9.97 

239.9 

-5.31 x lO- 3 

.3 

0 

6.2 

(b) 

4700 

! 




I 


; 9.99 

242.0 

j -5.46 X lO" 3 

i .1 

t 

! .8 

| 9.2 

(c) 

5500 

. 

! 

i 



! 

1 

9.97 

240.5 

,-5.00 x lO- 3 

1 

.3 

| .2 

0 

(d) 

6500 


i 

1 N 




; io.o 

240.4 

i -4.97 X lO" 3 

i 0 

.2 

.6 

10(a) 

3500 


i 

320 


; 

10.0 

240.0 

-4.84 X lO" 3 

0 

! 0 

3.2 

(b) 

4700 







9.94 

240.5 

-5.48 X lO" 3 

.6 

.2 

9.6 

(c) 

5500 





' 


9.97 

241.5 

-5.03 x lO" 3 

.3 

.6 

.6 

(d) 

6500 

V/ 

\ 

/ 

V/ 

31.1 

1634.6 

-4.55 X lO' 15 

211.0 

581.1 

100.0 



TABLE HI.- SUMMARY OF RESULTS FOR VARIED INITIAL VALUES OF T' 



Wavelength, 

A 

Initial conditions 

Final conditions 


Error, percentage 

Figure 

p', 

mb 

T . 

K 

r '> 

km-1 

p, mb 

T, K 

r, km-1 

El : - 10 x loo 
10 
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Figure 5.- Convergence plots for p’ = 18 mb, T’ = 240 K, and P = -5 X 10"^ km - l. 
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Figure 7.- Convergence plots for p’ = 10 mb, T' = 160 K, and I” = -5 X 10" 3 km" 1 . 
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Figure 8.- Convergence plots for p* = 10 mb, T' = 200 K, and F = -5 X 10"^ km -1 . 
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Figure 10.- Convergence plots for p* = 10 mb, T' = 320 K, and F = -5 X 10"^ km"*. 
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Figure 11.- Convergence plots for p’ = 10 mb, T’ = 240 K, and I” = +5 X 10“^ km - l. 
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Figure 12.- Convergence plots for p' = 10 mb, T' = 240 K, and T’ = 0. 



r; 

P,' 

t; 

km' 1 

mb 

K 

-8 X10- 3 

32 

500- 

r 

-7 

28 

r* 

450 

r 

-6 

24 

r 

400 — 


DESIRED 

VALUES 

O P 

□ T 

a r 


-5 

-4 





ITERATION NO. 

(c) X = 5500 A. 

Figure 12.- Continued. 


cn 

-a 



U1 

00 


T 

P, 

km' 1 

mb 

-32 X10- 3 

32 

-28 

28 

-24 

24 

-20 

20 

-16 

16 

-12 

12 

-8 

8 

-4 

4 

0 

0 


T, 

K 

500 
450 
400 
350 h 
300 - 


DESIRED 

VALUES 


O 

□ 

A 


P 

T 

r 


250 J— 

&- 


-rt 

A A 


A n o b~5 o'cTO ■cr cm-Q n-cro -e o- 

C ^ -o -Q - O O - Q g g O O O O o O O g-f O-g-i 0 - 0-0 0 Q O O O O O O 0 0 0 0 - 0 


200 


; o £ a □ a “ 

ba n o D ° a ° 


a A 


A A A AA A A 


150 fjA 


5 3 13 17 21 25 

ITERATION NO 

(d) X = 6500 A. 

Figure 12.- Concluded. 


j— A A A A A A A A A ft A 


1 nriA- -i l-i— i— L-i- 1 i I J- J I -I--J-I- i 1 i i i I i i i I i i i I i i i I i i i 

1 UUl C Q 1 O in? 0 1 rsr* r\ o oo 


29 33 37 41 


DESIRED 

VALUES 


r; 

p; 

X 

km" 1 

mb 

K 

-80 X10- 3 

32 

500- 

-70 

28 

450- 

-60 

24 

400 - 

-50 

20 

350- 

-40 

16 

300 - 


O P 

□ T 

a r 


-30 


-20 


-10 


0 


ct-B-B o-o-o -a- 0-o-o b-d-q -q- o-o-a o-o-o -e- q-o-b b-d-q -s- o-o-a b-o-q -e- q-o-b b-b-o 
QoOQOoe - OOOOOOOOQOOOOOOOOOO O O OOQOqOoOOO O o - o 



ITERRT ION NO. 


cn 

CO 


(a) \ = 3500 A. 

Figure 13.- Convergence plots for p' = 10 mb, T' = 240 K, and P = -8 X 10“^ km _ l. 
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Figure 13.- Continued. 
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Figure 13.- Concluded. 







DESIRED 

VALUES 


r ' p; t; 

km" 1 mb K 

-80 X 10- 3 32 500 - 


-70 

28 

450 - 

-60 

24 

400 - 

-50 

20 

350 - 


O 

□ 

A 


p 

T 

r 


-40 

-30 

-20 

-10 

0 


16 

12 

8 


300 - 


0 


250 - 


CJ-CI-B d-Q-O .0. Q-a-d 0-Q-a -D- B-O-B O-d-D -B- D-d -a B-D-B -o- B-O-B B-d-Q -& Q-d-Q B-D-d 
C POOOOOOOOOOO O ' O O O o ~Q oooooooooooooooo 0-0 o o o o-o 


200 
1 50 
100 


_aa AA aaaaaaaaaaaaa 


AAAAAAAAAA 


" 4" ' a" ' .V ' ‘ iV 1 ' g ' l' ' ' l ' ' ' as ' ' ' A ' 1 ' a 1 / ' 1 ' A 

ITERATION NO- 


05 

CO 


(a) X = 3500 A. 

Figure 14.- Convergence plots for p' = 10 mb, T' = 240 K, T' = -12 X 10“^ km"*. 
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Figure 14.- Continued. 
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Figure 14.- Continued. 
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Figure 15.- Convergence plots for p' = 10 mb, T' = 240 K, and r* = -50 X 10”^ km~^. 
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